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ABSTRACT 

Aims. We aim at directly detecting the presence of optically thin circumstellar dust emission within the terrestrial planetary zone 
around main sequence stars known to harbour cold debris discs. The present study focuses on a sample of six bright A- and early 
F-type stars. 

Methods. High-precision interferometric observations have been obtained in the near-infrared K band with the FLUOR instrument 
installed on the CHARA Array. The measured squared visibilities are compared to the expected visibility of the stellar photospheres 
based on theoretical photospheric models taking into account rotational distortion. We search for potential visibility reduction at short 
baselines, a direct piece of evidence for resolved circumstellar emission. 

Results. Our observations bring to light the presence of resolved circumstellar emission around one of the six target stars (f Aql) at 
the 5cr level. The morphology of the emission source cannot be directly constrained because of the sparse spatial frequency sampling 
of our interferometric data. Using complementary adaptive optics observations and radial velocity measurements, we find that the 
presence of a low-mass companion is a likely origin for the excess emission. The potential companion is characterised by a K-band 
contrast of four magnitudes. It has a most probable mass of about 0.6M o and is expected to orbit between about 5.5 AU and 8 AU from 
its host star assuming a purely circular orbit. Nevertheless, by adjusting a physical debris disc model to the observed Spectral Energy 
Distribution of the f Aql system, we also show that the presence of hot dust within 10 AU from f Aql, producing a total thermal 
emission equal to 1.69 ± 0.31% of the photospheric flux in the K band, is another viable explanation for the observed near-infrared 
excess. Our re-interpretation of archival near- to far-infrared photometric measurements shows however that cold dust is not present 
around ( Aql at the sensitivity limit of the IRS and MIPS instruments onboard Spitzer, and urges us to remove f Aql from the category 
of bona fide debris disc stars. 

Conclusions. The hot debris disc around Vega (Absil et al. 2006) currently remains our only secure resolved detection within the 
context of this survey, with six genuine early-type debris disc stars observed so far. Further observations will be needed to assess 
whether f Aql also belongs to this hot debris disc category. 

Key words. Stars: fundamental parameters - Circumstellar matter - Binaries: close - Methods: observational - Techniques: interfer- 
ometric 



1. Introduction 

Debris discs are optically thin, gas-poor dust discs around main 
sequence (MS) stars. The presence of circumstellar dust around 
stars with ages above ~10 Myr is attributed to populations of 
planetesimals that were neither used to make up planets nor 
ejected from the system by the time the nebular gas was dis- 
persed dMann et al.ll2006l) . These leftovers produce dust by mu- 
tual collisions and comet-type activity. Being continuously re- 
plenished by small bodies, the disc can then persist over much of 
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the star's lifetime. Due to its large total cross-section area, dust 
is much easier to observe than planets, not to speak of planetesi- 
mals. On the other hand, distributions of dust respond to the pres- 
ence of planetary perturbers, reflect distributions of the parent 
bodies and bear important memory of the planetary formation 
process in the past. Hence debris discs can be used as sensitive 
tracers of planets, as well as small body populations, and should 
reflect evolutionary stages of planetary systems. This explains 
the substantial effort invested in the observation and modelling 
of debris discs over the last two decades. 

Early-type stars, and A-type stars in particular, have been 
the most successfully studied targets in the context of debris 
disc studies so far. This is due to a great extent to their intrin- 
sic brightness, which efficiently lights up their debris discs up 
to large distances (Earth-like temperatures occur at about 5 AU 
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from such stars). Additionally, their MS lifetimes (~800 Myr) 
are long enough to encompass the main evolutionary stages of a 
typical planetary system. Far-infrared surveys carried out with 
space satellites have been particularly successful in detecting 
cold dust around nearby A-type stars. For instance, Spiter/MIPS 
observations of 160 A-type MS stars have shown that about 33% 
of such stars po ssess significant excess emission at 24 and 70 fim 
dSu et alj |2006). which represents a considerably higher excess 
rate than what has b een found for old solar analogs and M dwarfs 
dBrvden et alj|2006l) . 

Unlike these successful detections of cold debris in the far- 
infrared, the discovery of warm (~300K) or hot (~1000K) dust 
through mid- and near-infrared observations has been limited to 
a very small number of targets so far. Furthermore, the rare de- 
tections have generally been obtained for young MS stars show- 
ing strong mid-infr ared silicate features, su ch as ft Pic (~ 12 Myr , 
iPantin et all fl997|). HP 145263 (~ 8Myr, iHondaet all [2004), 
77 Tel (~ 12 M yr Ichenet all 120061) or HD 172555 (-12 Myr, 
IChen et alj[2006h . At such ages, planetary systems are still ex- 
pected to be in the process of forming planets, especially in the 
inner part of the disc (< 10 AU) where rocky planets may take 
up to one hundred Myr to accrete most of the small bodies on 
nearb y orbits and reach their final mass (see lKenvon & Bro mley 
2006,' for a solar-mass star). The observed large mid-infrared ex- 
cesses are therefore expected to be related to the end of the planet 
building phase rather than to the dust produced by the collisional 
grinding of an evolved planetary system, similar to the zodiacal 
dust in our solar system. For "mature" A-type stars, the lack of 
near- to mid-infrared excess has generally been associated with 
a dearth of warm dust grains, a suggestion generally confirmed 
by the presence of inner h oles in resolved images ar ound fa- 
mous objects such as a PsA dKalas et al.l2005l) or a Lyr dSu et all 
120051) . 

Incidentally, the inner part of debris discs is a very interest- 
ing region to study, as it directly probes the location where plan- 
ets are supposed to have formed and evolved — conversely, the 
outer part of the disc (> 20 AU), similarly to the Kuiper belt in 
our solar system, only bears a memory of the outermost massive 
planets through gra vitational interactio ns such as mean motion 
resonances (see e.g. lReche et alj |2008). The characterisation of 
inner debris discs is therefore crucial for an understanding of 
the formation, evolution and dynamics of planetary systems (in- 
cluding our own solar system), as well as to set the scene for the 
emergence of life on rocky planets. 

Reaching a better sensitivity to the inner part of debris 
discs has thus been an important challenge during the past 
years, and the detection of small amounts of hot dust around 
mature early-type stars has finally been enabled by the ad- 
vent of high-precision near-in frared stellar in t erfero metry. Using 
the VLTI/VINCI instrument, foi Folco et all d2004l) derived up- 
per limits to the Zf-band dust emission around five stars. The 
first robust resolved detection of hot d ust was obtained f or the 
bright AOV-type star a Lyr (Vega) by lAbsil et ail 



d2006l) with 

CHARA/FLUOR on an optimised set of interferometric base- 
lines, showing the presence of extended emission accounting for 
1.29 + 0.19% of the photo spheric flux in the K band. This result 
was already suggested by ICiardi et"aT] d200ll) . although with a 
large uncertainty on the flux ratio. 

This recent detection has raised a number of questions re- 
garding the nature and origin of inner dust grains. In particu- 
lar, a scenario involving the presence of star-grazing comets, in- 
jected into the inner planetary system by dynamical perturba- 
tions caused by migrating planets, has been proposed, inspired 
by the Late Heavy Bombardment (LHB) that happened in the 



early solar system (Gomes et al. 2005). A similar scenario is also 
proposed bv lWvatt"e t al. ( 2007) to explain the presence of warm 
dust (~300K) detected by the Spitzer Space Telescope around 
a few solar-type stars. We have thus decided to initiate a near- 
infrared interferometric survey of nearby debris disc stars to as- 
sess the occurrence of hot excesses around MS star s. The first 
paper of this series dDi Folco et alj|2007l hereafter iPaper J) fo- 
cused on two solar-type stars, showing the presence of hot dust 
around rCet. In this paper, we discuss early-type stars, which 
hold a privileged position in our survey because their brightness 
makes them well suited for an interferometric study. 

2. Methodology and stellar sample 

The principle of debris disc detection by stellar interferometry is 
based on the fact that the stellar photosphere and its surround- 
ing dust disc have different spatial scales. For an A-type star at a 
distance of 20 pc, the angular diameter of the photosphere is typ- 
ically about 1 mas, while the circumstellar disc extends beyond 
the sublimation radius of dust grains, typically located around 10 
to 20 mas for black body grains sublimating at r su b - 1500 K. In 
the infrared K band, the debris disc is therefore generally fully 
resolved at short baselines (10 - 20 m), while the photosphere 
is only resolved at long baselines (~200m). One can take ad- 
vantage of this fact to isolate the contribution of circumstellar 
dust by performing visibility measurements at short baselines, 
where the stellar photosphere is almost unresolved. The presence 
of resolved circumstellar emission then shows up as a deficit of 
squared visibility with respect to th e expec ted visibility of the 
bare stellar photosphere, as shown in Pa per j 



f(i)Ml-2teETO, 



(1) 



where 'V 2 and ( V\ are respectively the squared visibility of the 
star-disc system and of the bare stellar photosphere, b is the inter- 
ferometer baseline length, and ecsE is the flux ratio between the 
integrated circumstellar emission within the field-of-view and 
the stellar photospheric emission. This equation is valid only for 
short baselines and for ecsE <K 1, as expected for an optically 
and geometrically thin circumstellar disc. 

Pue to the expected faintness of the circumstellar emission 
(fcsE ^ 1 %X detection needs both a high accuracy on the mea- 
sured 'V 2 and on the estimation of the squared visibility 'V 2 of 
the bare photosphere at short baselines. While the former has al- 
ready been demonstrated in th e framework of sing le-mode near- 
infrared interferometers (e.g., Kervella et al. 2003), the latter is 
ensured to a large extent by the fact that the stellar photosphere 
is almost unresolved at short baselines, so that we can tolerate a 
certain level of imprecision on the knowledge of the stellar an- 
gular diameter without jeopardising the detection. For instance, 
a stellar photosphere with an angular diameter of 1 .0 ± 0. 1 mas 
produces a squared visibility of 0.995 ± 0.001 at a baseline of 
20 m in the K band. Therefore, in the following study, diameter 
measurements at long interferometric baselines are not required 
to derive accurate p hotospheric models, and empirical surface- 
brightness relations dKervella et all 2004) can instead be used to 
predict stellar angular diameters with a sufficient accuracy, tak- 
ing into account rotational distortion. 

To start our survey of early-type debris disc stars, we have 
chosen six bright (K < 4) nearby (~20pc) dwarfs, with spec- 
tral types ranging from AO V to F2 V. The main characteristics 
of our six target stars are su mmarised in Ta ble [1] together with 
the already-discussed a Lyr (Absil et al. 2006). All of them have 
been classified as debris disc stars based on the measurement of 
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a mid- or far-infrared excess emission above the expected photo- 
spheric level, using infrared satellites such as IRAS, ISO or the 
Spitzer Space Telescope. 

3. Observations and data reduction 

Interferometric observations were obtained in the infrared K 
band (1.94 - 2.34 fim) wit h FLUOR, the Fiber Linked Unit 
for Optical Recombination (ICoude du Foresto et al.l 120031) . us- 
ing the shortest baseline of the CH ARA Array formed by the S I 
and S2 telescopes (34 metres long, ten Brummelaar et al. 2005). 
Observations took place during Spring 2006, between April 30th 
and May 11th. 

The FLUOR field-of-view, limited by the use of single-mode 
fibres, has a Gaussian shape resulting from the overlap inte- 
gral of the incomin g turbulent w ave fronts with the fundamental 
mode of the fiber (lGuvonll2002l) . Its actual size depends on the 
atmospheric turbulence conditions. The estimation of the mean 
Fried parameter (ro) provided by the CHARA array during our 
observations (7 cm in the visible) allows us to derive a mean 
value of 0'.'8 for the full width at half maximum of the field- 
of-view. This parameter is sensitive to the actual atmospheric 
conditions for each individual target star (it typically ranges be- 
tween (X'7 and 1'.'2), but the only associated effect is to widen the 
search region for circumstellar emission. 

The FLUOR Data Reduction Software (DRS 
Coude du Foresto et al jl997t IKervella et al.l2004tlMerand et all 
2006) was used to extract the raw squared modulus of the co- 
herence factor between the two independent apertures. The 
extraction of the squared visibilities from the fringe packets 
recorded in the time domain is based on the integration of the 
squared fringe peak obtained by a Fourie r transform of the 
fringe packet (Cou de du Foresto et al.|[l 997). or equivalently by 
a wavelet analysis dKerveila et al.l 120041) . The interferometric 
transfer function of the instrument was estimated by observing 
calibrator stars before and after each observation of a scientific 
target. All calibrator stars (listed in Table |2]i were chosen from 
two catalogues developed for this specific purpose (Bord tTet al.l 
l2002t IMerand et al.ll2.005h . Calibrators chosen in this study are 
late G or K giants, whereas our target stars have spectral types 
between AO and F2. Since the visibility estimator implemented 
in the FLUOR DRS depends on the actual spectrum of the 
target star, an appropriate correction must be applied to our 
data, otherwise our squa red visibilities would be bia sed at 
a level of about 0.3% (ICoude du Foresto et al l 1 19971) . This 
cor rection can be based ei t her on the shape factors discussed 
by ICoude du Foresto et al. I (119971) or on a wide band model 
for estimating the calibrator's visibilities and interpretin g the 
data (see e.g. iKervella et al.l l2003; Auf denberg et al.l l2006). The 
latter method was chosen for this work, and all the calculations 
presented here therefore take into account a full model of the 
FLUOR instrument, including the spectral bandwidth effects. 

4. Estimating stellar angular diameters 

A reliable estimation of the photospheric angular diameter is an 
important pre-requisite for the detection of circumstellar emis- 
sion at short baselines. In this section, we use surface-brightness 
relations to estimate the mean photospheric diameter of our tar- 
get stars. We then apply a standard model of a rotating pho- 
tosphere in hydrostatic equilibrium to estimate the rotation- 
induced distortion. Our approach is finally validated by com- 
parison with long-baseline interferometric measurements where 
available. 



Table 2. Ca librators chosen from the catalogues of iBorde et al.l 
(120021) and IMerand et ail d2005l) . listed with spectral type, K 
magnitude, uniform disc (UD) angular diameter in K band with 
lcr error bar, and associated scientific target(s). 



Identifier 


Sp. type 


m K 


0w 


(mas) 


Targets 


46 Boo 


K2III 


2.66 


1.433 


±0.019 


ttCrB 


9QB 


G9III 


3.07 


1.101 


±0.014 


aCrB 


QY Set 


K8 Illb 


1.53 


2.680 


± 0.029 


aCrB 


HD 92095 


K3 III 


2.20 


1.527 


± 0.020 


/JUMa 


HD 100615 


K0III 


3.19 


1.027 


±0.014 


/JUMa 


44UMa 


K3 III 


2.04 


1.990 


± 0.023 


/JUMa 


X Vir 


K2III 


2.02 


1.960 


± 0.023 


r\ Crv 


HD 108522 


K4III 


3.46 


1.171 


±0.016 


r\ Crv 


HD 111500 


K4III 


3.05 


1.395 


± 0.020 


j] Crv 


HD 157617 


K1III 


2.71 


1.305 


±0.018 


y Oph 


710ph 


G8 III 


2.43 


1.500 


± 0.020 


7 Oph, 












a Boo, £ Aql 


HD 162113 


K0III 


3.62 


0.904 


±0.012 


y Oph 


HD 162468 


Kl III-IV 


3.17 


1.123 


±0.015 


y Oph 


HD 133392 


G8 III 


3.26 


1.052 


±0.014 


cr Boo 


HD 132304 


K3 III 


3.48 


1.042 


±0.014 


cr Boo 


HD 126597 


K2III 


3.46 


1.004 


±0.014 


cr Boo, y Oph 


HD 176527 


K2III 


2.04 


1.721 


± 0.024 


i Aql 


HD 175743 


K1III 


3.35 


1.106 


±0.015 


( Aql 


f Aql 


G9.5 Illb 


2.37 


1.620 


±0.021 


^ Aql 


/" Aql 


K3 Illb 


1.76 


2.240 


± 0.023 


f Aql 



4.1. Diameter estimation from surface-brightness relations 

Accurate angular diameter estimations for MS stars can be ob- 
tained by applying empiric surface-brig htness relations (see e.g., 
IKervella et alJ2004tlDi B enedetto 2005), using as input the mea- 
sured magnitudes in two standard photometric bands. These re- 
lations have been c alibrated by interferometric measurements 
(IKervella et al] l2004). and are valid for spectral types from AO 
to M2. The smallest intrinsic dispersion (cr < 1%) is obtained by 
combining K- or L-band magnitudes with visible magnitudes (B 
or V). Here, we choose to use the /f-band magnitudes, which are 
generally known with good accuracy, and com bine them with 
the V -band magnitudes measured by Hipparcos ( Perrvma n et al.l 
1997). K magnitudes have been preferred to L magnitudes, be- 
cause the potential presence of hot dust around the target stars is 
expected to affect the K band at a smaller level than the L band — 
the contrast between a ~9000K photosphere and ~1000K sur- 
rounding dust decreases with increasing wavelength. In prac- 
tice, we do not expect the Zf-band magnitudes to be affected 
by dust emission by an amount larger than 5%, as for higher 
contributions the signature of hot dust would have already been 
detected. Our interferometric measurements will confirm this as- 
sumption a posteriori. The (V, V- K) relation, where the V mag- 
nitude is used to calibrate the reference magnitude, has been cho- 
sen in the present study because of its robustness to dust emis- 
sion and because the accuracy of the V magnitudes (0.01 maj], 
iPerrvman et al.l 119971) is better than that of the K magnitudes. 
The final error bar on the estimated limb-darkened angular di- 
ameter #ld is therefore generally smaller using the (V, V - K) 
relation than the (K, V - K) relation. The mathe matical expres- 
sion o f the (V, V-K) relation is (see Table 4 from lKervella et ail 
12004 : 

log 9ld(V, V -K) = 0.2753(V-Zf) + 0.5 175 -0.2V. (2) 



1 This estimated error bar is conser vative (or even pessimistic ) for the 
bright stars that we are studying here ( Ivan Leeuwen et afl l997). 
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Table 1. Fundamental parameters and estimated angular diameters for the seven selected targets (including a Lyr). A mean limb- 
darkened diameter is first computed using Eq. [2] — this actually represents the geometric mean of the minor and major axes of the 
elliptical photosphere. The apparent photospheric oblateness is then estimated by means of Eq.[5] from which we finally deduce the 
angular diameters for the apparent minor and major axes of the photosphere. The Icr errors are given in superscript. 



Name 


HD 


Type 


Dist. 


Mass 


v sin i 


m v 


m K 


Mean 6 LD 


Appar. 


Minor #ld 


Major $ld 


Refs. 








(pc) 


(M ) 


(kms" 1 ) 






(mas) 


oblat. 


(mas) 


(mas) 




ySUMa 


95418 


Al V 


24.3 


2.28 


46 


2.35 UU1 


238 u.U6 


1.095 u " 2i 


1.007 




1.099 l)uii 


1, 4, 6, 7 


n Crv 


109085 


F2V 


18.2 


1.44 


92 


4.30 001 


3.54 005 


0.736 0015 


1.022 


0.728 0015 


O.7440015 


1, 4, 6, 8 


cr Boo 


128167 


F2V 


15.5 


1.28 


15 


4.47 001 


3.49° 02 


0.783 0016 


1.001 


0.783 0016 


0.783 0016 


1, 4, 6, 9 


aCrB 


139006 


AOV 


22.9 


2.58 


138 


2 ; 23 o.oi 


2.20° 05 


1.202 0024 


1.059 


1.1 68° 024 


1 237 0024 


2, 4, 6, 10 


7 Oph 


161868 


AOV 


29.1 


2.18 


210 


3 75 01 


3.67 005 


0.616 0012 


1.107 


0.585 0012 


0^648 0012 


1, 4, 6, 11 


a Lyr 


172167 


AOV 


7.8 


2.30 


22 


oio3 001 


0.00 002 


3 312 067 


1.002 


3 309 0067 


3.3 15° 067 


3, 5, 6, 12 


<TAql 


177724 


AOV 


25.5 


2.37 


317 


2.99 001 


2.90° 02 


o'880 0018 


1.307 


; 770 o.oi 8 


1.006° 018 


1,4, 6, 13 



Refer enc es. — Masses fro m (l)[Allende Pri eto & LamberJ 19991. (2)lTomkin & Popped 19861, (3)|Auldenberg et al. 2006; v sin i f rom (4)lRover et al 



2002, (5) 



1996, (9) 



Hill e t"ai]2004 V ma gnit udes from (6l | Perrvman et al.ll 19971: K mag nitud es from (7)lNeugebauer & L eighton 1997, (8) Sy lvester et al. 
Blackwellet all 19791 (lO) ISneden et all 19781 (ll) ISelbv et al.lll988l , U2) lAumann & Probsj|199lL (13) lLeggett et alj|1986l 



It has a dispersion of 1 % and shows no significant non linearity. 
The angular diameter estimations for our seven targets are given 
in the ninth column of Table Q] ("Mean #ld")- 

The associa t ed lim b darkening coefficients uk can be found 
in IClaret etaD (1 19951) . assu ming a linear l i mb-da rkening law 
hip) — 1 - «i(l - fj) as in iKervella et al.l (|2004|) . where p is 
the cosine of the angle between the perpendicular to the surface 
of the star and the line of sight. The closest tabulated models to 
the physical properties of the stars were chosen. It must be noted 
that the rapid rotation of some of our targets may lead to sub- 
stantial gravity darkening, so that the tabulated limb darkening 
coefficients could be underestimated. Fortunately, considering a 
photosphere 1.0 mas in diameter observed with a 30 m baseline, 
even a threefold increase of the limb darkening coefficient would 
produce a variation smaller than 0.1% on the squared visibility 
(see Eq. [6] in Sect. 0. The case of a Lyr, which is slightly re- 
solved on a 30-m baseline, must however be considered with 
additional care (see Sect. 14. 3t . 

Estimating the angular diameter of a CrB is a little more 
complicated than for the other six targets, as it is actually a bi- 
nary system of Algol type. The primary and secondary are es- 
timated to be AO V and G5 V stars respectively, and the mag- 
nitude difference in V b and is estimated to be AMy = 4.89 
(Tom kin & Popper! 119861) . Based on a V-band magnitude of 
2.22 + 0.01 for the system (fPerrvman et al.|[T997l) . we adopt a 
V-band magnitude of 2.23 + 0.01 for the prima ry. Then, using 
the in trinsic colours of MS stars in the infrared (Bes sel & Brett! 
1988), we compute the magnitude difference between the pri- 
mary and secondary in the K band (AM^ = 3.25), and derive the 
expected /f-band photometry of the primary from the system's 
photometry in a similar way. The original K magnitu de for the 
binary system, K = 2.14 + 0.04 dSneden et al.lll978l) . gives a 
corrected magnitude of 2.20 ± 0.05, where the error bar has been 
slightly revised to account for possible errors in the estimation of 
the spectral types and intrinsic colours of the two components. 



surface of the ellipse to be equal to the angular surface derived 
from surface-brightness relations in the previous section: 



R a Rb - R„ 



8, n d 



(3) 



where R m is the radius of the equivalent spherical photosphere 
(i.e., the geometric mean of the major and minor radii) and d the 
distance to the target star. Under the assumptions of hydrostatic 
equilibrium, uniform rotation, and of a point mass gravitational 
potential, the apparent rotational velocity v sin i can be related to 
the major and minor radii R a and R/, with the following equation 
dTassoullll978l) : 



v sin 1 



2GM 



1 



1 



Rb Ra 



(4) 



We define the apparent oblateness (p) of the photosphere as the 
ratio of the major and minor radii (p = R a /Rt,), and solve the 
above equation for p as a function of R m to obtain: 



p =* 



(vsiru) R„ 
4GM 



1 + 



(v sin i) 2 R„ 
4GM 



(5) 



Using the limb-darkened angular radii, masses, and v sin i listed 
in Table[T] we derive the apparent oblateness of the photospheres 
of our seven target stars. They range between 1.001 and 1.307 
(see Table [U. From these estimations of oblateness, we can de- 
rive the expected angular sizes of the major and minor axes. 
Since the orientations of the photospheres are unknown, the 
differences between the major and minor angular radii directly 
give the additional uncertainty on the stellar diameters along the 
CHARA S 1-S2 baseline due to stellar oblateness. The additional 
error is typically a few percent of the mean stellar diameter, and 
reaches 14% for the fastest rotator. 



4.2. Effect of rapid rotation 

With projected rotational velocities up to 317kms _1 (Table Q]), 
we expect the photospheres of our target stars to deviate signif- 
icantly from spherical symmetry due to large centrifugal forces. 
We will assume in the following study that the apparent pho- 
tospheric shape can be modelled by an ellipse with semi-major 
and semi-minor axes noted R a and Rh, and we force the angular 



4.3. Long-baseline interferometric measurements 

Since long-baseline infrared interferometric measurements are 
generally not available for our target stars, we will try to val- 
idate our approach by comparison with interferometric mea- 
surements of similar MS stars. For the first step of our mod- 
elling approach, i.e., surface-brightness relation s, such a vali- 
dation has already been thoroughly addressed by IKervella et al.l 
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Table 3. Application of surface-brightness relations and of our standard rotating photosphere model to four rapidly rotating stars 
that have been directly measured with infrared stellar interferometers. The agreement between our estimated mean limb-darkened 
stellar diameter ("$ld model") and oblateness ("p model") with the measured stellar parameters is very satisfactory. The largest 
discrepancy between model and observations is fo und for a Aql and is less than 3% — it should be noted that the interferometric 
measurements reported on this particular target by iMonnier et alj d2007l) and[Domiciano de Souz a et al. I J2QQ5h disagree at the 3cr 
level (4% of LD ). 



Name 


Type 


v sin i 


Mass 






P 


P 


Ref. 






(kms" 1 ) 


(M ) 


model 


meas. 


model 


meas. 




a Leo 


B7V 


317 


3.4 


1.47 


1.47 


1.323 


1.320 


(1) 


a Lyr 


AOV 


21.9 


2.3 


3.31 


3.31 


1.0015 


1.0014 


(2,3) 


a Cep 


A7 IV-V 


283 


2.0 


1.51 


1.54 


1.296 


1.298 


(4) 


a Aql 


A7V 


240 


1.8 


3.31 


3.39 


1.170 


1.164 


(5) 



Refere nces: (1) iMcAlister etaf] J2005l) . (2) lAufdenberg et al.1 d2006f) . (3) iPeterson et all d2006bl) . (4) Ivan Belle et al] J2006I) . (5) IMonnier et ail 
d2007h . 



(120041) . Now, for the second step, i.e., the elliptic photospheric 
model, we can test the validity of Eq. [5] on four rapidly rotat- 
ing stars with similar masses to our target stars, whose apparent 
oblateness has been derived from interferometric measurements 
(see Table [3]). The agreement between the modelled apparent 
oblateness ("p model") and the observed oblateness ("p meas.") 
is very good. 

To improve the accuracy of our angular diameter estima- 
tions, actual long baseline interferometric measurements would 
however be valuable. In the case of a Lyr, such measurements 
are available, and have been used to produce a realistic model of 

I Ilk 

the stellar photosphere in the K band ( Aufdenberg et al. 2006). 
Using the K-bwA limb profile of lAufdenberg et all (120061) . we 
have derived an accurate limb-darkened model for the stellar 
photosphere, giving the following parameters: angular diame- 
ter #ld = 3.305 + O.OlOmas and linear limb-darkening coeffi- 
cient uk - 0.3610 The angular diameter is sligh tly smaller than 
the ph ysical diameter (3.33 mas) derived by lAufdenberg et all 
(2006) because the outermost layers of the stellar atmosphere, 
from 0.9927?, to R, (where the pressure is less than 0.01 
dynes/cm 2 ), does not contribute significantly to the flux emis- 
sion. The angular diameter obtained from surface-brightness re- 
lations (3.312 + 0.067 mas, see Table [TJ is in good agreement 
with this measurement, although it must be noted that the associ- 
ated error bar is 7 times larger than that of the directly measured 
diameter. Nonetheless, we will see in Sect. [5] that the accuracy 
on the modelled angular diameters displayed in Table [TJ is gen- 
erally sufficient to allow the detection of circumstellar emission 
at a level as low as 0.1%. The case of a Lyr must be consid- 
ered separately, because its 3.305 mas photosphere is partly re- 
solved on a 30 m baseline in K band. The error on the modelled 
photospheric diameter would then induce a large enough uncer- 
tainty on the stellar 'V 2 (about 0.6%) to compromise the detec- 
tion of circumstellar emission. Moreover, due to gravity darken- 
ing, the actual limb-da rkening coefficien t is much larger than the 
tabulated one (0.20 in lClaret et al.lll995l) . Therefore, in the rest 
of this study, we will us e for a Lyr the photospheric model of 
Aufdenb erg et al.l ((2006) described here above. 

Another target in our sample has been observed with long- 
baseline interferometry, this t ime in the visible rang e: £ Aql was 
used as a reference star by IPeterson et al.l (l2006al) to calibrate 
interferometric observations of Altair with the Navy Prototype 



2 The limb-darkening coefficient varies significantly across the K 
band. We have computed here a broad-band coefficient integrated on 
the FLUOR pass band. 



Optical Interferometer (NPOI). Using archival NPOI data, the 
authors derived a preliminary model for the £ Aql photosphere, 
giving a minimum angular diameter of 0.815 ± 0.005 mas. Even 
though this value is well within our angular diameter range 
(see Table [TJ, it is 2.5cr from the minimum angular diameter 
that we derive from our rapidly rotating photospheric model 
(0.770 ± 0.018 mas). This discrepancy could partly be due to the 
wavelength difference between the observations (V band) and 
the model (K band). 

We also note th at <x Boo has been observed as a calibrator by 
Bode n et"ai] d2005) at the Palomar Testbed Interferometer (PTI) 
and at the NPOI. The authors adopted an angular diameter of 
0.77 ± 0.04 mas based on SED modelling of broadband photom- 
etry, a value which is consistent with our estimation. They did 
not mention any significant deviation from that estimated value 
based on their measurements. 

Finally, the eclipsing bin ary nature of a CrB provid es an ad- 
ditional check of our model: iTomkin & Popper! (1 1 986b derive a 
limb-darkened radius of 3.04 + 0.03 R e for the primary, which 
translates into an angular diameter of 1.235 + 0.012 mas. This 
value is most probably representative of the major axis of the 
photosphere as the rotation axis of the primary and the orbital 
axis of the binary system are expected to be aligned. The agree- 
ment with our major angular diameter (1 .237±0. 024 mas) is very 
good. 

5. Searching for circumstellar emission 

To assess the presence of circumstellar emission around a tar- 
get star, we compare the interferometric 'V 2 measurements with 
the expected squared visibility of the stellar photosphere, using 
the oblate photospheric model presented in the previous section. 
Assuming a linear limb darkening law, the monochro matic stel- 
lar visibility is given by dHanburv Brown et alJI 19741) : 

2 ll-u A u A Y 2 1 .Ji(x) u A ^h /2 (x)\ 2 

^•* = (— -jj l (1 "^-r + ^r-^j ' (6) 

where x = n6bb/A, where b is the interferometer baseline length 
projected on the sky plane and 6b is the limb-darkened diame- 
ter along the direction of the baseline. The wide-band visibilities 
are then computed by integrating the above expression over the 
FLUOR bandp ass, using the stellar spectrum as a weighting fac- 
tor (see Pa per J) . 

The unknown orientation of the stellar rotation axis is taken 
into account by considering two extreme cases: the first case 
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("large diameter") assumes that the major axis of the elliptic 
stellar photosphere is aligned with the mean orientation of the 
interferometric baseline projected onto the sky plane for all ob- 
servations, while the second case ("small diameter") assumes the 
minor axis to be aligned with the mean baseline orientation. We 
further add the 1 <x uncertainty on the mean limb-darkened diam- 
eter (see TableQ} to the large diameter case, while we subtract it 
from the small diameter case. The first case leads to the small- 
est visibilities, and the second case to the largest visibilities. In 
that way, we produce a "lcr box" for the squared visibility of 
the photosphere at any projected baseline orientation (or equiv- 
alently, any hour angle of observation on the S1-S2 baseline), 
which is represented by the (thick) blue lines in Fig.Q] The com- 
bination of good precision on the stellar model and a small base- 
line length — which leaves the photosphere almost unresolved — 
results in a very small region of possible values for the squared 
visibilities of the naked stars. The small curvature of the lcr box 
is due to the fact that the projected baseline length is maximum 
at an hour angle of h and decreases before and after meridian 
crossing. 

5.1. Data quality and fitting procedure 

The calibrated 'V 2 have been overlaid on the plots of the ex- 
pected photospheric visibilities in Fig.Q] The quality of the data, 
which can be estimated by the size of the error bars and the dis- 
persion of the data points, is not uniform for all targets, and de- 
pends mostly on the stellar magnitude and on the seeing condi- 
tions during the observations. The low SNR obtained on indi- 
vidual measurements for the faintest targets (rj Crv, <x Boo and 
y Oph) leads to the largest error bars (typically ranging between 
2% and 4% for these targets, depending on the particular condi- 
tions of the observations). The error bars for the brightest targets 
typically range between 1% and 2%. The statistical dispersion of 
the data points shows a healthy behaviour for all targets, as will 
be confirmed by the values of the reduced^ 2 obtained when fit- 
ting the data with a simple star+disc model. The relatively large 
scatter in the a CrB data is due to poor seeing conditions (esti- 
mated ro - 4.5 cm) on the night this target was observedQ 

The calibrated 'V 2 data have been fitted with a model con- 
sisting of a limb-darkened stellar photosphere surrounded by a 
uniform circumstellar emission filling the whole field-of-view 
of the instrument. To account for stellar oblateness in our fit, we 
have defined for all baseline orientations a mean photospheric 
diameter as the arithmetic mean of the angular diameters in the 
"large diameter" and "small diameter" cases (which both depend 
on the orientation of the projected baseline). The only parame- 
ter to be fitted is then the flux ratio ecsE between the integrated 
circumstellar emission and the photosphere. The results of the 
fitting procedure are displayed as insets in Fig. Q] together with 
the reduced^ 2 . The first term in the error budget accounts for the 
statistical dispersion of the data, while the second term is related 
to the uncertainty on the photospheric model, which includes the 
uncertainties on both the surface-brightness model and the ori- 
entation of the stellar rotation axis. 

5.2. Results 

Among the seven early-type stars that have been surveyed so far, 
two stars {( Aql and a Lyr) show evidence for the presence of 
circumstellar emission: the fit to their interferometric measure- 

3 The mean seeing for the other five targets of the May 2006 observ- 
ing run is ro - 8 cm. 



ments falls significantly below the expected squared visibility 
of the photosphere (see Fig.Q] and Fig. [2j. The five remaining 
stars in our sample show no evidence for circumstellar emission: 
the estimated flux ratio ecsE between the integrated circumstellar 
emission and the photosphere ranges between -0.3<x and +0.9<x, 
where <x represents the square sum of the statistical and system- 
atic errors bars displayed in Fig. Q] It must be noted that our 
non-detections have a statistically-speaking healthy behaviour, 
with a mean A"-band excess amounting to 0.3<x and a standard 
deviation of 0.6<x. This gives us good confidence in our positive 
detections, which are obviously not related to a systematic effect 
since our stellar sample is quite homogeneous in spectral types 
and magnitudes. The range of values for the reduced x 1 of the 
fits is satisfactory for all cases. Let us discuss the most interest- 
ing targets individually. 



£ Aql. Our analysis of the new CHARA/FLUOR data brings to 
light the presence of circumstellar emission around ( Aql (see 
Fig.Q}. The quality of the fit for ( Aql is quite good (j^ 2 = 0.97), 
and the final value for the £"-band flux ratio is 1.69% + 0.31%, 
which represents a 5.5<x detection. This result does not strongly 
depend on the actual limb darkening coefficient uk of the ( Aql 
photosphere: when using extreme values of 0.0 and 0.5 for uk, 
the final value for the jf-ba nd excess varies o nly from 1.68% 
to 1.72%. As discussed by Absi l et all (120061) . this estimation 
of the /f-band flux ratio does not depend on the morphology of 
the circumstellar emission to the first order, as long as it has an 
extension larger than the angular resolution of the interferometer 
(about 7 mas in the K band for a 34-m baseline). 



cCrB. Let us assess whether the eclipsing binary nature of the 
object could significantly affect our measured visibilities. To re- 
duce as much as possible the effect of the companion on the mea- 
sured visibilitie s, we have used the well-characterised ephemeris 
of the system dSchmitt) 1 1 998i) to schedule the observations of 
a CrB during an eclipse. As discussed in Sect. 14.11 the expected 
/T-band magnitude difference between a CrB A and its compan- 
ion is AK = 3.25, which is equivalent to a contrast of about 5%. 
The companion could thus potentially be at the origin of the fluc- 
tuations of about 10% observed in the calibrated squared visibil- 
ities. While the secondary eclipse would have been ideal to sup- 
press the effect of the companion, time and weather constraints 
have forced us to observe the system during primary eclipse 
(when the companion passes in front of the primary). The a CrB 
system was observed on April 30th from 8:10UT to 12:04 UT, 
while the primary eclipse took place between 6:56 UT (start of 
ingress) and 21:19 UT (end of egress). Totality of the eclipse ac- 
tually started at 9:36 UT, so that part of the observations were 
obtained during ingress. Assuming that the companion is a stan- 
dard G5 V star (see Sect. 14. It , and using for the primary the linear 
radius of 2.94 R G deduced from the surface-brightness relation- 
ships, the diameter ratio between the two components is equal to 
0.31. With the estimated A'-band contrast, the expected variation 
of squared visibility for the binary system on the S1-S2 baseline 
between the start of ingress and the totality is only about 0.6%. 
This effect is clearly not sufficient to explain the dispersion of the 
data points, which is more likely related to the poor atmospheric 
conditions during the night of April 30th. Due to our optimised 
observation timing, the presence of the companion does not bias 
the search for a circumstellar excess at a level larger than about 
0.2% on the disc/star contrast. 
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Fig. 1. Squared visibilities with error bars as a function of hour angle, recorded for our six targets on the CHARA S1-S2 baseline. 
The blue region represents the lcr box defined by our photospheric models. We have fitted to our data a model of a limb-darkened 
photosphere surrounded by a uniform circumstellar emission. The best-fit model is represented by the black solid line, while the 
dotted lines represent the lcr statistical uncertainty on the best-fit model. Circumstellar emission is detected when the blue region 
lies significantly above the uppermost dotted line. The inferred contrast between the uniform circumstellar emission and the star is 
given as inset for each star ("disc/star contrast"), together with its two error bars (the first one is related to the statistical dispersion 
of the calibrated data set and the second one to the uncertainty on the photospheric model) as well as the reduced^- 2 of the fit. 



77Crv. The absence of significant K-band excess around this (Wva tt et"ai] 120051). and subsequent modelling an LHB-like 
target comes rather as a surprise, since previous mid-infrared event dWvatt et al.ll2007l) . However, we note that (i) the latest 
observations have suggested the presence of warm dust model of the 77 Crv circumstellar disc predicts a contrast larger 
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Fig. 2. Same as Fig.Q]in the case of a Lyr. The previously pub- 
lished data have been re-reduced with the latest version of the 
FLUOR Data Reduction Software to produce this plot, using the 
latest limb-darkened photospheric model of lAufdenberg et al] 
(2006) as a reference. 



than 1:10 s in the K band between the disc and the photosphere 
(Smith et al. 2008) and (ii) the small number of data points that 
we have obtained on this target does not allow a high sensitivity 
to be reached (only 1:62 at 3<x). Therefore, our non-detection is 
not incompatible with previous measurements and models, and 
shows that the presence of warm dust is not necessarily accom- 
panied by large amounts of hot dust. 

a Lyr. The data collected in 2005 on a Lyr have been re- 
reduced with the latest version of the FLUOR DRS, and using 
a more realistic model f or the stellar photosphere based on the 
Aufdenbe rg et al.l d2006l) model (a simple uniform disc model 
was used in the original study bv lAbsil et alj|2006l) . The new fit, 
displayed separately in Fig. [2] shows a good agreement with the 
previously published /T-band excess (1.29% ± 0.19%). Only the 
error bar on the disc/star contrast has significantly changed, due 
to the slightly larger error bars obtained on individual data points 
with the latest version of the FLUOR DRS. The contribution of 
the uncertainty on the photospheric model is negligible (0.06%) 
thanks to the accurate stellar model derived from long-baselines 
measurements. The agreement between the result presented here 
and the original disc/star contrast estimation can be seen as a 
welcome consistency check, since a new version of the FLUOR 
DRS was used to reduce the data and a more realistic stellar 
model was adopted. This also shows that the details of the pho- 
tospheric model do not matter at the short baselines considered 
here. 



6. Possible nature of the circumstellar emission 
around f Aql 

The nature of the circumstella r emission detected around a Lyr 
has already been discussed bv lAbsil et alj d2006l) . showing that 
the presence of hot circumstellar dust is the most probable expla- 
nation. Here, we focus on the new detection (£ Aql), and extend 
our analysis to other potential contributors to circumstellar emis- 
sion, including stellar winds and mass loss events. The following 
discussion is based on our interferometric observations, comple- 
mented only by archival infrared photometric measurements. In 



Sect. 13 we will include additional data from other instruments 
to further constrain the nature of the circumstellar emission. In 
this way, we want to make clear what type of information can 
actually be derived from a limited amount of squared visibilities 
obtained on a short two-telescope baseline. 

The only strong constraints that we derive on the circum- 
stellar emission from our interferometric measurements is that 
it represents 1.69% of the stellar flux in K bancQ and that it is 
located within the FLUOR field-of-view, which has a FWHM 
of 078. Due to the quasi-Gaussian beam of the FLUOR single- 
mode fibres, circumstellar emission located at the edge of the 
field (e.g., dust ring or point-like source) needs to be signif- 
icantly brighter than 1.69% to produce the detected visibility 
deficit. Taking the transmission at the centre of the field as a 
reference, only half of the flux emitted at 074 (i.e., at a projected 
distance of 9 AU) makes it through the fibre's transmission pat- 
tern. At 30 AU, the flux transmission decreases to 10%, follow- 
ing the theoretical Gaussian profile. 

6.1. The debris disc scenario 

The first scenario that we wish to investigate is the presence of 
hot circumstellar dust within the field-of-view of the FLUOR in- 
strument, which would be at the origin of the measured visibility 
deficit. Although the origin and survival of hot dust grains close 
to an early-type MS star is a matter of debate, a plausible de- 
bris disc model, reproducing the interferometric measurements 
as well as archival photometric measureme nts in the near- and 
mid-infrared has been proposed for a Lyr (lAbsil et al.ll2006l) . In 
that paper, we also showed that the actual morphology of the 
circumstellar dust disc does not have a significant influence on 
the visibility deficit measured by CHARA/FLUOR: only the flux 
ratio between the integrated disc emission and the stellar photo- 
sphere actually matters. 

In this section, we investigate whether a debris disc model 
similar to that proposed for a Lyr would be suitable to repro- 
duce the available measurements of £ Aql. The purpose is not to 
propose a "best-fit" set of physical parameters for the dust disc, 
but only to check whether there exists at least a debris disc model 
consistent with the near- and mid-infrared observations of f Aql. 
The question of the origin and physical properties of such dust 
grains will be discussed in a forthcoming paper. 

To perform our simulations, we have complemented our 
direct measurement of the /T-band excess flux with archival 
spectro-photometric measurements at near- and mid-infrared 
wavelengths (see Table @). In addition to the measurements 
found in the literature, we have used an IRS spectrum found 
in the Spitzer archives , which we have re duced with the c2d 
pipeline developed by Lah uis et al.l ((2006). The spectrum has 
been binned into a few equivalent broad-band photometric mea- 
surements for the sake of SED modelling (see Table|4]i. Because 
all these fluxes are valid for the whole star/disc system, they 
have been converted into an estimated excess flux by subtract- 
ing the photospheric flux. The p hotospheric emission m odel is 
chosen from the NextGen grid (Hauschild t et al.|[l999l) . using 
stellar para meters as close a s possib le to the pa r amete rs es- 
timated by iGrav et al.l (l2003l) and lAdelman et ail d2002l) . i.e., 
r eff = 9200 K, logg = 4.00 and M/H = 0.0. The NextGen 
model is scaled to match the observed V = 2.988, which re- 
sults in an esti mated luminosity of 39L B , i n agreement with the 
estimation of iMalagnini & Morossil (I1990I) . We have estimated 
the error on the NextGen photospheric infrared fluxes by com- 



The effective wavelength of the FLUOR observation is 2. 132 jim. 
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Table 4. Available constraints on th e near- and mid-infra r ed pho - 
tome try of 4" Aql . Refe rence s: (1) iMorel & Magnenatl d 19781): 
(2) iLeggett et al.l dl986l) : (3) lEgan et all d2003b : (4) IChen et all 
(120051) . The archival Spitzer-IRS spectrum has been binned into 
a few photometric data points. Photometric data have been com- 
pared to our NextGen photospheric model of ^"Aql to derive 
fractional excesses, taking into account the typical uncertainty 
of 2% on the photospheric flux in the infrared (this uncertainty 
has been added to the estimated errors on the measured fluxes to 
compute the errors on the excesses). 



Wavelength 


F 

1 meas 


(Jy) 


Excess 


Ref. 


1.2 yum 


112.4 


±4.1 


-4.5 


± 5.5% 


(1) 


1.24 yum 


114.5 


+ 2.1 


2.0 


t 3.9% 


(2) 


2.13 ytim 






1.69 


±0.31% 


CHARA/FLUOR 


4.35 yum 


13.3 d 


- 1.7 


-2.3 


± 14.7% 


(3) 


5.5 /im 


9.18 ± 


0.42 


5.3 


± 6.8% 


IRS spectrum 


7.5 /im 


4.85 ± 


0.22 


1.7 


± 6.6% 


IRS spectrum 


8.28 yum 


4.34 ± 


0.18 


9.5 


± 6.5% 


(3) 


12.13pm 


2.15 ± 


0.13 


14.8 


± 8.9% 


(3) 


12.5 yum 


1.73 ± 


0.10 


-1.8 


± 7.6% 


IRS spectrum 


14.65 fim 


1.44 ± 


0.10 


12.1 


± 9.8% 


(3) 


18yum 


0.868 ± 


0.051 


2.0 


± 8.0% 


IRS spectrum 


24yum 


0.475 ± 


0.048 


-0.7 


± 12.0% 


(4) 


32yum 


0.270 ± 


0.025 


0.8 ± 


11.5% 


IRS spectrum 



paring neighbouring models on the NextGen grid corresponding 
to the typical range in T e ff and \ogg estimations found in the 
literature. Assuming effective temperatures comprised between 
9200 K and 9400 K, and surface gravities between 3.5 and 4.0, 
we have derived a typical error of 2% on the estimated fluxes 
between 1 fim and 32 yum. 

Our debris disc model is an adaptat ion of the "best- fit" model 
that we deriv ed in the case of or Lyr dAbsil et al.l l2006). based on 
the model of Augereau et al. ( 1999J- It assumes a size distribu- 
tion dn(a) cc a~ 3 - da with limiting grain sizes a mm =0.1 yum 
and a max = 1500 /im, a surface density power-law S(r) oc r -4 , 
and chemical composition of 50% amorphous carbon and 50% 
glassy olivine. The inner radius of the dust distribution has been 
recomputed in a similar way as in the case of a Lyr, and is found 
to be equal to 0.14 AU. This corresponds to the sublimation ra- 
dius for the grains larger than about 1 /zm, while 0. 1 yum grains 
sublimate at ~ 0.45 AU in the model. All grains within their sub- 
limation distance are self-consistently eliminated from the sim- 
ulations, leading to regions around the inner disc edge depleted 
in submicronic grains. The only free parameter is the total mass 
of the disc, which has been adjusted to the photometric and in- 
terferometric data in Fig[3] With a reduced^ 2 of 0.54, our model 
does a good job in reproducing the global SED of f Aql from 1 
to 32 /mi, and demonstrates that the presence of hot dust, with a 
total dust mass of only 5.5 x 10~ 8 M ffl and a fractional luminos- 
ity of about 10~ 3 , is a viable explanation to the observed /f-band 
excess. It must be noted that, as in the case of a Lyr, most of the 
grains in this model are submicronic and located close to their 
sublimation radius. Possible scenarios to explain the presen ce of 
such grains have already been proposed dAbsil et al.ll2006l) and 
will be further discussed in a forthcoming paper. 

The hot debris disc model proposed for f Aql can further 
be used to derive upper limits on the mass of hot dust around 
the stars in our sample for which no detection was reported. 
This estimation is based on the fact that none of the surveyed 
stars shows a significant photometric infrared excess in the 1- 
5 /im region, and on the fact that their spectral types are close 
enough for the same debris disc model to apply to all of them. 



Noting that the final error bar on the disc/star flux ratio is typi- 
cally 0.5% for these stars (except for /3 UMa, which has a signif- 
icantly lower error bar of 0.16%), we derive a typical 3cr upper 
limit of 5 x 10~ 8 M ffi for the hot dust content within the FLUOR 
field-of-view of 078 (i.e., 18 AU at the mean distance of our sam- 
ple). Scaling this mass to the smaller error bar of ySUMa gives 
a 3<t upper limit of 1.5 x 10~ 8 M© for this particular system. 
We note that our observations are only sensitive to hot dust, and 
that larger amounts of colder dust could be present, especially 
beyond a few AUs where equilibrium temperatures are signifi- 
cantly smaller than 1000 K (see e.g. the discussion of ^Crv in 
Sect.|521). 



6.2. Stellar wind and mass loss 

With their relatively high masses and high rotational velocities, 
one may think that our target stars could be the subject of signif- 
icant mass loss or winds that could, as in the Be phenomenon, 
lead to a substantial near-infrared emission in addition to the 
photospheric level. However, winds of A-type stars are actually 
suspected to be very weak, as the mechanisms which are effi- 
cient in both later type stars and earlier ones are not expected to 
operate on A-type stars: radiative pressure is very small (com- 
pared to B-type stars), and there is no deep con vection givin g 
rise to activity phenomenon as in late-type stars dBabelll 1995h . 
Convection is actua lly expected to start for e ffective tempera- 
tures below 6500 K dLamers & Cassinellil 19991) . i.e., for spectral 
types later than about F6 V, so that A-type stars are not expected 
to have chromospheres and coronae. This fact seems to be cor- 
roborated by the lack of observational evidence for winds ema- 
nating from A-type stars dLanz & Catalal[l992l lLandstreet et al.l 
1998), as we ll as by the absence of X- ray emission from most 
A-type stars dSchroder & Schmittll2007l) . 

Nonetheless, we note that the presence of a stellar wind has 
been detected around Sirius (Al V) through the observation of 
an absorption f eature in the blue wings of the Mg II and H I 
resonance lines dBertin et al.ll 19951 ). The inferred mass-loss rate 
ranges between 2xl0~ 13 and 1.5xlO _12 M yr _1 , and has not been 
shown to produce any significant n ear-inf r ared e mission. Based 
on radiatively driven wind models, Babel (1995) predicts much 
lower mass loss rates (< 10 _I6 M o yr _1 ), with winds consisting 
of only metals. Additionally, the infrared excess that could be 
associated to the winds from hot s tars is expected to originat e 
mostly from a region R < 1.5/?. dLamers & Cassinellil fl999). 
which does not match our observations because such a compact 
emission would not produce a significant visibility reduction at 
the short baselines used here. 

Intriguingly, shell lines of Tin have been detected by 
lAbt et al.l d 19971) in the ultraviolet spectra of approximately 
one-quarter of the most rapidly rotating normal A-type dwarfs 
(v sin i > 200 km s _1 ). In that paper, the authors suggest that this 
"hot disc" phenomenon might be related to sporadic mass-loss 
events and does not seem to correlate to the prese nce of cold 
dust in an outer debris disc. However, more recently, Eb3 {2004) 
and lAbtl J2008) have argued that these hot discs are actually ac- 
creted from the interstellar medium. The argument is based on 
the fact that such discs rarely occur around stars within the heart 
of the Local Interstellar Bubble (LIB). The author proposes a 
model in which stars accrete discs in dense interstellar regions, 
but loose them in regions of low interstellar density, such as the 
LIB, where the stellar winds exceed the accretion rate. The spec- 
tral lines associated with the hot disc disappear on timescales of 
roughly decades. Since our target stars are all located well within 
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Fig. 3. Top panel: A possible fit of our debris disk model to the photometric and interferometric constraints of Table|4] The dashed 
line represents the contribution of thermal emission to the overall disc emission, while the solid line includes the scattered light 
contribution. Photometric measurements are represented by squares with error bars, and the related excesses by diamonds. The 
CHARA/FLUOR interferometric measurement is represented by a filled circle with error bar. The whole IRS spectrum is represented 
by a grey box with a width equal to its error bar, and was converted into a few regularly sp aced (in log scale) broad-band photometric 
measurements for the fitting procedure. The MIPS photometric measurement at 70jum dChen et al.ll2005l) is also represented but 
was not included in the fit. Botto m panel: Same figure in the case of a low-mass companion. The spectrum of the companion is 
represented by a NextGen model (Hauschil dt et al.ll99 9) where the stellar parameters (r e ff = 3800 K, logg = 4.5) have been chosen 
as close as possible to our best-fit model (see Sect. 16.31 ). The resolution of the NextGen spectrum of the AO-type host star has been 
degraded by a factor 20 in both panels for the sake of clarity. 



the LIB (<K lOOpc), this phenomenon does not seem appropriate 
to reproduce our observations. 

All these elements indicate that a stellar mass-loss origin 
to the observed near-infrared excess is highly unlikely in the 
case of A-type dwarfs. Furthermore, we have examined archival 
HARPS data on £ Aql and found no sign of Ha emission, which 
accompanies strong stellar winds produced by hot stars. The stel- 
lar wind scenario will thus be discarded in the rest of the discus- 
sion. 

6.3. Point-like source 

Because of our sparse sampling of spatial frequencies, we cannot 
determine the actual morphology of the excess emission source, 
and a point-like source located within the FLUOR field-of-view 



(either a bound companion or a background source) could be at 
the origin of the observed visibility deficit. To assess the effect of 
a faint companion on the measured visibilities, two cases must 
be considered, depending whether the two fringe packets are su- 
perimposed or not. Taking into account the coherence length of 
the FLUOR fringes (about 25 ftm), the minimum on-sky angular 
separation for two separated fringe packets is 150 mas along the 
baseline direction, i.e., about 3.5 AU at the distance of £ Aql. 
For larger separations, the off-axis point-like source contributes 
essentially as an incoherent emission and produces a squared 
visibility deficit A'V 2 = 2/, where / is the flux ratio between 
the primary and the secondary. A point-like source within the 
FLUOR field-of-view, located at a distance larger than 150 mas 
from f Aql in the projected direction of the interferometric base- 
line and producing 1.69% of the photospheric flux of £Aql, 
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Fig. 4. Map of the reduced x 2 distance between a binary star 
model and our CHARA/FLUOR squared visibilities for £ Aql, 
as a function of the position of the faint companion. The binary 
flux ratio has been adjusted separately at each given position to 
minimise the x 2 - Dark regions indicates binary solutions com- 
patible with the CHARA/FLUOR data. 



arated fringe packets can also be used as a good approximation 
for superposed fringe packets. 

Figure|4]is also a nice illustration of the wide bandwidth ef- 
fect in optical interferometry: for position angles parallel to the 
mean baseline orientation (i.e., PA = 170°), one can see that 
the ripples in the x 2 map due to the visibility modulation pro- 
duced by the binary star progressively wash out as the compan- 
ion moves away from the primary star. This effect can be inter- 
preted as the progressive separation of the two fringe packets. 
The reduced x 2 values obtained at a position where the fringe 
packets are essentially separated (i.e., separation of 150 mas and 
position angle close to 170°) make a nice connection with the 
case of large separations described earlier in this section. One 
can also note that for position angles perpendicular to the mean 
baseline azimuth, the reduced^ 2 does not change significantly as 
a function of distance, because only the projection of the angular 
separation along the baseline direction actually matters. 

Finally, we have checked that we can find on the NextGen 
stellar atmosphere grid a low-mass stellar companion compat- 
ible in terms of flux both with our CHARA/FLUOR measure- 
ment and with the archival photometric data listed in Table [4] In 
the bottom panel of Fig. [3] we show that an MOV-type star with 
T e ff = 3800 K and \ogg = 4.5 nicely fits the whole data set. A 
low-mass companion located at an angular distance larger than 
3 mas from the central star is therefore a plausible explanation to 
the observed visibility deficit at short baselines. 



would therefore be a possible explanation to the measured visi- 
bility deficit. 

Let us now investigate the case of a close companion (< 
150 mas), for which the two fringe packets would be superposed. 
In that case, the companion produces a modulation in the fringe 
visibility depending on its position and on the geometry of the 
array, which changes as the Earth rotates. The squared visibili- 
ties are therefore expected to be modulated as a function of the 
hour angle of the observation. To determine which combination 
of angular separation and position angle would be compatible 
with the interferometric observations, we compute the squared 
visibilities for a binary system at various distances and position 
angles, assuming that the off-axis source did not move during the 
1 1 nights of our May 2006 observing run0 For each given posi- 
tion, the contrast between the two components is then adjusted 
to minimise the x 2 distance between the synthetic and measured 
squared visibilities. In this way, we produce ax 2 map, taking into 
account the wide bandwidth effect of the FLUOR instrument. 

The x 2 map, represented in polar coordinates in Fig. |4] 
indicates that many different locations of the off-axis point- 
like source may reproduce the squared visibilities measured by 
CHARA/FLUOR. In particular, for any given angular separa- 
tion larger than 3 mas, we can always find a position angle that 
makes the simulated data compatible with the observed visibili- 
ties (reduced ;f 2 < 1). We also note that the optimum binary flux 
ratio at the positions where the reduced x 2 is smaller than 2 is 
in excellent agreement with the "disc/star" contrast computed in 
Sect. 16.11 the mean flux ratio at such locations is 1.65%, with 
a standard deviation of only 0.22%. This further confirms that 
the morphology of the circumstellar emission does not have a 
significant influence on the inferred flux ratio, and indicates that 
the simplified computation of the binary flux ratio assuming sep- 



7. Discussion 

Our interferometric observations are not sufficient to determine 
the morphology of the circumstellar emission. In this section, we 
investigate whether complementary data obtained with various 
types of instruments could further constrain the nature of the 
circumstellar excess. 



7. 1 . Further constraints on the debris disc scenario 



3 This hypothesis is valid for background sources and for bound com- 
panions with semi-major axes larger than 0.7 AU (which does not pre- 
vent the companion from being at a projected angular separation as 
small as a few mas from the target star). 



The id entification of f Aql as a debris disc star by IChen et all 
(2006) was based on the measurement of significant excess 
emission on top of the expected photospheric flux at 12 and 
24yum. However, Fig. [3] strongly contradicts this statement, and 
shows the absence of cold circumstellar dust at the sensitivity 
limit of the Spitzer I RS and MIPS instru ments. The discrepancy 
with the analysis of IChen et alj ([2006) is related to the effec- 
tive temperature of 1 1912 K that they have used to compute the 
photospheric emission of £Aql, which does not match the es- 
timated temperatures fo und in the literature, ranging between 
9 190 and 9680 K (see e.g.lGrav et al.l207)llAdelman et all2002t 
Malagnini & Morossl| [l990|). Only the polar tem perature could 
reach such a high value (Peterso n et al.1 1200 6a). while £Aql is 
view essentially equator-on. The very good match between our 
photospheric emission model and the IRS spectrum gives us a 
high confidence in our flux estimation, and clearly shows that 
no excess emission is detected in the 24-70/mi region . This also 
explai ns why a black body disc model did not allow IChen et al.l 
d2005l) to consistently reproduce the 24 and 70 yum excesses that 
they had derived. The fact that the MIPS 70 pm measurement 
falls about 3cr below our photospheric model (see Fig. [3]) could 
be either due to a poor flux calibration on this target or to an 
underestimated error bar. 

Despite the absence of observable amounts of cold dust, 
the presence of hot circumstellar dust in the close neighbour- 
hood of 83 Myr old A-type star could still be explained in 
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the context of terrestrial-plan et formation models. In particu- 
lar, iKenvo n & Bromley (2006) describe oligarchic and chaotic 
growths in the terrestrial region around a solar-type star, and 
show that the planet formation timescale is 10-100Myr (see 
also IChambersI l200l and references therein). The onset of 
chaotic growth depends on the surface density, and typically 
ranges between 0. 1 Myr and a few Myr. During the subsequent 
chaotic phase, the dynamical perturbations between oligarchs 
may lead to enhanced dust production. This potential origin for 
the dust observed a round ( Aql is also supported by the study of 
Riek e et ail (|2005), who suggest that the dust detected towards 
relatively young A-type stars may be related to the planet accre- 
tion end game. According to lCham bers (2001), the accretion end 
game is predicted to be dominated by a small number of colli- 
sions between large bodies (e.g., the early Earth and the impactor 
that caused format ion of the Ea r th's M oon) and may stretch 
over 100-200 Myr. iKokubo etaf] d2006l) also show that the gi- 
ant impact stage (i.e., the accretion timescale) around a Solar- 
type star last s for about 100 Myr. Howev er, it must be noted that, 
according to IKenvon & B romley (2004}), the major part of the 
warm infrared excess may disappear within 1-10 Myr during 
terrestrial planet formation around Sun-like star. Furthermore, 
IKenvon & Bromlevl (120051) show that, in the case of a 3-20 AU 
region around A-type stars, 1000 km objects are formed within 
10 Myr, the dust luminosity peaks at 1 Myr and is already at a 
very low level at 100 Myr. 

In conclusion, even though the presence of hot dust in the 
close environment of (Aql cannot be ruled out, the likeliness 
of this scenario can be questioned due to the absence of a large 
reservoir of cold dust in the outer planetary system. 

7.2. Further constraints on the binary scenario 

To get a better grasp at the possible nature of the potential point- 
like source around ( Aql, we compute the flux of the hypothet- 
ical companion reproducing the observed visibility deficit as a 
function of its distance, taking into account the Gaussian off- 
axis transmission of the single-mode fibre. For the sake of clar- 
ity, the apparent binary flux ratio seen through the Gaussian 
transmission pattern is approximated by the simple expression 
(/ = A r V 2 /2), which strictly applies only to wide binaries 
(> 150 mas), even though we have shown in Sect. 16. 3 1 that this 
relationship is a good approximation to the best-fit solutions for 
close binaries. The derived flux ratio is then converted into the 
mass of the hidden companion in Fig.|5](solid line) as a function 
of linear distan ce to the cent r al star , using the mass-luminosity 

relationships of iBaraffe et all (Il998l). 

As already discussed by I Absil et al.l (120061) and iPaperll the 
presence of a background source producing such a flux within 
the 078 FLUOR field-of-view is highly unlikely. Therefore, we 
investigate here the possible presence of a bound companion 
around £ Aql. Although ( Aql is known to be part of a wide mul- 
tiple system (ADS 12026), its detected companions are no closer 
than 5" so that they cannot affect our visibility measurements. 
This star is classified as a spectroscopic binary in the WDS cat- 
alogue. However, this classification is most probably wrong as it 
is based on a measurement from the late 19th century, which has 
not been confirmed since then. The detection of X-ray emission 
towards ( Aql with ROS AT is a possible sign for th e presence 
of a low-mass companion ( Schroder & Schmitt 2007). However, 
due to the large PSF of the ROSAT satellite (18"), the visual 
companion that was located at 6'.'5 from ( Aql in 1934 according 
to the WDS catalogue could be at the origin of the observed X- 
rays. To our knowledge, no other sign of the presence of a close 
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Fig. 5. Constraints on the mass and distance of a hypothetic com- 
panion around £ Aql. The solid line, labelled FLUOR, represents 
the solutions that are compatible with our interferometric obser- 
vations. The adaptive optics images that we have obtained with 
PUEO constrain the companion to be located below the dashed 
line (PUEO 3<x detection limit). The probability that the com- 
panion defined by the FLUOR curve fits our HARPS data is rep- 
resented by three colour levels along the FLUOR curve. The cor- 
respondance between the mass of the companion and its TT-band 
contrast with respect to the primary is indicated by grey hori- 
zontal lines, based on the evolutionary models of IBaraffe et al.l 
(I1998I) . 

companion to ( Aql has been reported in the literature. However, 
close and faint companions are difficult to detect. 

In order to improve our knowledge of the close environment 
of £ Aql, deep adaptive optics images have been obtained in 
the K band with the PUEO instrument at the Canada-France- 
Hawaii Telescope. Under good seeing conditions (about 076), 
no companion was detected at the sensitivity limit of the instru- 
ment (AK = 4 at 073, and improving for larger separations). The 
3<x sensitivity to point-like companions has been computed as a 
function of angular separation and has been converted into the 
maximum mass of a hidden companion in Fig. [5] (dashed line) 
as a functio n of linear distance, using the mass-luminosity rela- 
tionships of lBaraffe et al.ldl998l) . All solutions with companions 
orbiting further than about 8 AU can be discarded by the PUEO 
observations, but closer companions are not constrained. 

We have therefore used radial velocity data obtained with 
HARPS to further constrain the possible presence of a low-mass 
companion around ( Aql. We have tried to fit simple orbital so- 
lutions to the HARPS data collected between 2004 and 2008 
(represented in Hig. [6j, assuming a 90° inclination for the binary 
system and a zero eccentricity. Each point along the FLUOR 
curve in Fig.|5]defines a single couple "period - maximum radial 
velocity", which is computed with classical orbital dynamics re- 
lationships and defines the shape of the sinusoidal radial velocity 
curve. For each couple "period - velocity" on the FLUOR curve, 
we have computed the^ 2 between the observed HARPS data and 
a series of sinusoidal models with various time and radial veloc- 
ity offsets. The lowest^- 2 , associated to the best-fit offsets in time 
and radial velocities, has then been converted into a probability, 
which represents the likeliness that the binary star model actu- 
ally reproduces the HARPS data. Probabilities are represented 
by three colour levels in Fig. [5] and show that two solutions at 
orbital radii a = 1.2 AU and a = 2.5 AU are most probable (i.e., 
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Fig. 6. Radial velocities measured on £ Aql with HARPS during 
the past 4 years. The two most probable circular orbits compati- 
ble with both the HARPS data and the FLUOR constraints have 
been overlaid as solid and dotted curves. 



would not have been detected at the 3<x level by Hipparcos. This 
corresponds to a semi-major axis of about 5.5 AU. All solutions 
with semi-major axes between 0.2 AU and 5.5 AU would thus 
most likely have been detected by Hipparcos, and the two most 
probable solutions provided by the fit to the HARPS data seem 
incompatible with the astrometric constraints. However, we note 
that there remains a range of plausible solutions at semi-major 
axes between 5.5 AU and 8 AU, with stellar masses ranging be- 
tween 0.6 and 0.65 M Q . We also note that the use of eccentric 
orbits to fit the HARPS data could lead to significantly different 
orbital solutions, which may improve the compatibility with the 
Hipparcos constraints. 

In conclusion, the presence of a bound low-mass companion 
can neither be confirmed nor rejected based on currently avail- 
able data. We plan to obtain further observations of £ Aql at high 
angular resolution with single-pupil telescopes to give a final 
answer to this question. In particular, the new Sparse Aperture 
Masking mode of the NACO adaptive optics instrument at the 
VLT would be perfectly suited to reach this goal. Interferometric 
closure phase measurements on longer baselines could also help 
unveiling the nature of this infrared excess. 



47 mas and 98 mas respectively). The potential companions have 
respective masses of 0.56 and 0.57 M© and complete their orbits 
in 280 and 842 days. These two most probable solutions are re- 
spectively represente d by solid and dotted curves in Fig. [6] The 
closest model on the Baraff e~et alj dl998l) evolutionary grid at 
80Myr for these two potential companions is a 0.57 M object 
with an effective temperature T e s = 3779 K and a surface grav- 
ity logg = 4.68. Its absolute magnitude is Mk = 5.22, giving 
K = 7.26 at the distance of ( Aql, and thus a AK = 4.36 with the 
primary. The associated V-band magnitude is 1 1 .2, so that the 
contrast with the primary is AV = 8.2. 

Even though the HARPS data do not give a definitive an- 
swer regarding the orbital parameters of the putative compan- 
ion, we must recognise that the absence of a bound companion 
seems rather unlikely based on the HARPS data alone: when fit- 
ting the HARPS data with a constant velocity, the reduced y 1 
with 9 degrees of freedom is equal to 3.6, and the associated 
probability that the model represents the data is of only 0.02%. 
Nevertheless, this result depends essentially on two HARPS ra- 
dial velocity measurements with values around 5000 m/s in the 
middle of Fig. [6] and a better time sampling would be needed to 
confirm the detection of the putative companion. It must finally 
be noted that stellar pulsations are not expected to produce such 
radial velocity variations in the case of £ Aql, a main-sequence 
AO star not classified as a 6 Scuti variable. 

A last constraint that we can put on the potential low-mass 
companion is that it must comply with the astrometric stability 
of ( Aql: Hipparcos observed this star on about 20 different oc- 
casions between March 1990 and March 1993 (total time span 
of 1000 days), and f ound no sign of binarit y at an astrometric ac- 
curacy of 0.72 mas (Perrv man et al.lll997l) . Taking an estimated 
mass of 0.57 M for the companion, this basically rules out 
any circular binary solution with a semi-major axis larger than 
0.2 AU at the 3<x level. However, if the orbital period is signif- 
icantly larger than the time span of the Hipparcos observations, 
the presence of a low-mass companion could have been missed, 
all the more that the system is supposedly viewed edge-on so 
that the reflex motion of the primary star could have been con- 
fused with its proper motion. We estimate that binary solutions 
with orbital periods larger than about 3000 days, even though 
they are associated with astrometric signals larger than 50 mas, 



8. Conclusion 

In this paper, we have investigated the close neighbourhood 
(< 1") of six nearby A- and early F-type MS stars, in 
search for hot counterparts to the cold debris discs detected by 
mid- and far-infrared spectro-photometric space-based observa- 
tions. The high-accuracy squared visibilities collected with the 
CHARA/FLUOR interferometer, combined with semi-empirical 
models of stellar photospheres including rotational distortion, 
has allowed us to reach dynamic ranges ranging from 1:175 
to 1:625 at lcr around the target stars. At this level of pre- 
cision, the presence of a resolved K band emission has been 
identified around only £ Aql, with a estimated Zf-band excess of 
1.69+0.31%. This detection ad ds to ou r previous results on a Lyr 
( Abs il et all 120061) and rCet dPaper J) , giving an overall near- 
infrared excess detection rate of 3/9 for the MS stars surveyed 
so far, among which 2/7 are early-type stars. The healthy statis- 
tical behaviour of the five non-detections in the present sample 
and the confirmation of the excess emission around a Lyr with an 
improved photospheric model and a new version of the FLUOR 
Data Reduction Software demonstrate the robustness of our ap- 
proach for hot debris disc detection. 

Our near-infrared interferometric measurements are not sam- 
pling the Fourier frequency plane in a sufficiently dense manner 
to derive the morphology of the excess emission source. In par- 
ticular, both a point-like source and an extended circumstellar 
emission can reproduce our observations. While in the cases of 
or Lyr and rCet, the presence of a bound or unbound companion 
to the target stars within the small FLUOR field-of-view could 
be rejected with a high confidence, we cannot rule out the pres- 
ence of a low-mass companion in the close vicinity of £ Aql to 
reproduce the measured K-bwd excess. The combination of the 
astrometric stability of ^"Aql measured by Hipparcos, the vari- 
ability of the radial velocities measured with HARPS and the ab- 
sence of off-axis companion in PUEO observations restricts the 
parameter space of the high-contrast binary scenario to masses 
in the range 0.6 to 0.65 M Q and semi-major axes between 5.5 
and 8 AU (i.e., about 200 and 300 mas). The ^-band contrast be- 
tween the primary and its companion would then be AK ~ 4, 
making it one of the closest high-contrast companions resolved 
around MS stars so far. 
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Besides a low-mass companion, the presence of hot circum- 
stellar dust grains producing a significant thermal emission in 
the K band is another viable explanation of the observed excess 
emission. In parti cular, we show that the debris disc model that 
Absi l et aTl (1200 6) have proposed in the context of ahyv is con- 
sistent with both our /T-band detection and archival near- and 
mid-infrared spectro-photometric measurements. However, our 
re-interpretation of archival Spitzer/MIPS and Spitzer/IRS data 
cle arly shows that th e presence of an outer debris disc, suggested 
bv lChen etail (f2005). can be firmly ruled out at the sensitivity 
level of MIPS and IRS. In the absence of significant amounts of 
cold dust, the hot debris disc scenario is not favoured to explain 
our CHARA/FLUOR measurements. 

The statistics of the "hot debris disc" phenomenon presently 
remains poorly constrained: for early-type stars, a hot debris disc 
was found around only one star (a Lyr) out of six bona fide de- 
bris disc stars observed so far. The case of a Lyr could therefore 
be rather unusual. Our interferometric survey of debris disc stars 
will be extended to larger sample in the coming years to improve 
our statistics. The case of £ Aql must be considered separately, 
since we show that it is not surrounded by cold dust at the sen- 
sitivity level of the Spitzer instruments, and since a close com- 
panion is a likely explanation to the observed /f-band excess. 
This star will deserve special attention in the future, and further 
observations will be performed to determine the actual nature 
of the near-infrared excess emission that we have resolved. In 
particular, infrared aperture masking experiments on large tele- 
scopes have the potential to reveal the true nature of the observed 
excess. 
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